Molecular markers for the assessment of genetic variability in threatened plant species

La variabilità genetica è un importante strumento per lo studio e la conservazione della biodiversità in specie rare e minacciate di estinzione. Durante il mio dottorato mi sono quindi occupata di mettere a punto diverse metodologie molecolari al fine di valutare la diversità genetica in due specie rare della flora italiana che presentano problematiche diverse e specifiche. I marcatori arbitrari RAPD e i marcatori semi-arbitrari ISSR sono stati utilizzati per valutare la diversità genetica in Quercus crenata Lam. e per confermare l’ipotesi della sua origine ibridogena dalle due specie presunte parentali Quercus cerris L. e Quercus suber L., essendo Q. crenata presente in Italia settentrionale dove Q. suber è attualmente assente. I marcatori SSR o microsatelliti sono invece stati messi a punto su una specie a rischio di estinzione, endemica dell’Appennino Tosco-Emiliano, Primula apennina Widmer, applicando una metodologia specifica, basata sulla costruzione di una libreria genomica arricchita per l’isolamento di primer specifici. I marcatori RAPD e ISSR, utilizzati su un totale di 85 campioni, hanno mostrato alti livelli di diversità molecolare entro le specie studiate, eccetto per Q. suber le cui popolazioni rappresentano il margine orientale di distribuzione della specie, per questo più sottoposte ad impoverimento genetico. Oltre alla cluster analysis (UPGMA) e alla Analisi delle Componenti Principali effettuate per entrambi i marcatori, che confermano l’ipotesi dell’origine ibrida degli individui di Q. crenata diffusi in Italia Settentrionale, sono stati calcolati l’indice di ibridità basato sul maximum likelihood, che dimostra una introgressione asimmetrica di Q. crenata verso il parentale caratterizzato da superiorità demografica (Q. cerris) e il test di Mantel. Quest’ultimo ha permesso di confrontare i due marcatori RAPD e ISSR utilizzati ottenendo una bassa correlazione, a conferma del fatto che, amplificando tratti differenti del DNA nucleare, i dati non sono sovrapponibili, sebbene forniscano risultati analoghi. Per l’isolamento di loci microsatelliti ipervariabili ho utilizzato il protocolllo FIASCO (Fast isolation by AFLP of sequences containing repeats- Zane et al. 2002) che permette di costruire una libreria genomica arricchita partendo dal DNA estratto da P. apennina. Tale procedura ha previsto la digestione del DNA genomico per la produzione di una miscela di frammenti di DNA. Tramite ibridazione con opportune sonde sono stati isolati i frammenti contenenti i microsatelliti. Sequenziando i cloni ricombinanti, ho ottenuto sequenze contenenti repeats sulle cui regioni fiancheggianti sono stati costruiti 15 coppie di primer che potranno, in seguito, essere utilizzate per definire la quota di riproduzione clonale in P. apennina e per valutare la diversità genetica delle popolazioni che coprono l’areale di distribuzione della specie. Data la loro natura altamente variabile e la loro abbondanza nel DNA, gli SSR saranno, come i marcatori RAPD e gli ISSR, ugualmente validi per lo studio della variabilità genetica e per l’analisi di problematiche specifiche legate alle specie rare.

TNFRSF13B Genetic variability an anthropological - evolutionary approach to Biomedical Research

In the recent years TNFRSF13B coding variants have been implicated by clinical genetics studies in Common Variable Immunodeficiency (CVID), the most common clinically relevant primary immunodeficiency in individuals of European ancestry, but their functional effects in relation to the development of the disease have not been entirely established. To examine the potential contribution of such variants to CVID, the more comprehensive perspective of an evolutionary approach was applied in this study, underling the belief that evolutionary genetics methods can play a role in dissecting the origin, causes and diffusion of human diseases, representing a powerful tool also in human health research. For this purpose, TNFRSF13B coding region was sequenced in 451 healthy individuals belonging to 26 worldwide populations, in addition to 96 control, 77 CVID and 38 Selective IgA Deficiency (IgAD) individuals from Italy, leading to the first achievement of a global picture of TNFRSF13B nucleotide diversity and haplotype structure and making suggestion of its evolutionary history possible. A slow rate of evolution, within our species and when compared to the chimpanzee, low levels of genetic diversity geographical structure and the absence of recent population specific selective pressures were observed for the examined genomic region, suggesting that geographical distribution of its variability is more plausibly related to its involvement also in innate immunity rather than in adaptive immunity only. This, together with the extremely subtle disease/healthy samples differences observed, suggests that CVID might be more likely related to still unknown environmental and genetic factors, rather than to the nature of TNFRSF13B variants only.

Phylogeny and genetic diversity of Italian species of hares (genus Lepus)

In this study we have analysed the genetic variability in ca. 700 samples belonging to six species of genus Lepus using maternal and biparental molecular markers (mitochondrial DNA, microsatellites, Single Nucleotide Polimorphisms). We aimed to reconstruct the phylogenetic relationships of species of hares living in Europe, and assess the occurrence of hybridization between the European hare Lepus europaeus and the Italian hare Lepus corsicanus. Results showed a deep genetic differentiation and absence of hybridization between L. corsicanus and L. europaeus, confirming that they are distinct and distantly related biological species. In contrast, we showed small genetic distances and a close phylogenetic relationship between the Italian hare and Cantabrian hare L. castroviejoi, which suggest a deeper evaluation of their taxonomic status. Populations of L. corsicanus are geographically differentiated. In particular, the peninsular and Sicilian populations of Italian hares are sharply genetically distinct, which calls for avoiding any translocation between Italy and Sicily. Information on genetic variability and population structure is being used to implement the Italian Action Plan for L. corsicanus.

Genetic characterization, population history and evolutionary medicine perspective in two native south american populations: Yanesha and Wichi

Two Amerindian populations from the Peruvian Amazon (Yanesha) and from rural lowlands of the Argentinean Gran Chaco (Wichi) were analyzed. They represent two case study of the South American genetic variability. The Yanesha represent a model of population isolated for long-time in the Amazon rainforest, characterized by environmental and altitudinal stratifications. The Wichi represent a model of population living in an area recently colonized by European populations (the Criollos are the population of the admixed descendents), whose aim is to depict the native ancestral gene pool and the degree of admixture, in relation to the very high prevalence of Chagas disease. The methods used for the genotyping are common, concerning the Y chromosome markers (male lineage) and the mitochondrial markers (maternal lineage). The determination of the phylogeographic diagnostic polymorphisms was carried out by the classical techniques of PCR, restriction enzymes, sequencing and specific mini-sequencing. New method for the detection of the protozoa Trypanosoma cruzi was developed by means of the nested PCR. The main results show patterns of genetic stratification in Yanesha forest communities, referable to different migrations at different times, estimated by Bayesian analyses. In particular Yanesha were considered as a population of transition between the Amazon basin and the Andean Cordillera, evaluating the potential migration routes and the separation of clusters of community in relation to different genetic bio-ancestry. As the Wichi, the gene pool analyzed appears clearly differentiated by the admixed sympatric Criollos, due to strict social practices (deeply analyzed with the support of cultural anthropological tools) that have preserved the native identity at a diachronic level. A pattern of distribution of the seropositivity in relation to the different phylogenetic lineages (the adaptation in evolutionary terms) does not appear, neither Amerindian nor European, but in relation to environmental and living conditions of the two distinct subpopulations.

The genetic history of Italians: new-insights from uniparentally-inherited markers

Genetic differences among human groups can be ascribed both to the broad-scale extents of pre-historical and historical migrations and to the fine-scale impacts of socio-cultural and geographic heterogeneity. In this thesis, the genetic information provided by uniparental markers were exploited to address different aspects of the Italian population history, by combining macro- and micro-geographic investigations at different spatial and temporal scales. To firstly assess the overall Italian variability, Y-chromosome and mtDNA markers were deeply typed in ~900 individuals from continental Italy, Sicily and Sardinia. Sex-biased patterns and contrasting demographic histories were observed for males and females. Differential European and Mediterranean contributions were invoked to explain the paternal genetic sub-structure observed in peninsular Italy, compared to the homogeneous maternal genetic landscape. If Neolithic showed to be one principal determinant of the detected paternal structure, local insights into specific Italian regional contexts highlighted the importance of Post-Neolithic contributions. Among them, migrations from the Balkans (particularly Greece) during late Metal Ages, played a relevant role in the cultural and genetic transitions occurred in Sicily and Southern Italy. On a finer geographic and temporal perspective, the more recent layers of Italian genetic history and some aspects of the gene-culture interaction were assessed by exploring the genetic variability within two “marginal populations”: Arbereshe of Southern Italy and Partecipanza in Northern Italy. The Arbereshe are Albanian-speaking communities settled in Sicily and Calabria since the end of Middle Ages. Despite sharing common genetic and cultural backgrounds, these groups revealed diverging micro-evolutionary histories, implying different founding events and different patterns of cultural isolation and local admixture. Partecipanza is an idiosyncratic institution of Medieval origin aimed at sharing and devolving collective lands. This case-study exemplified that socio-economic stratification within the same population may induce sex-biased genetic structuring and the maintenance of otherwise hidden historical genetic traces.

Innovative analytical methods for Central Nervous System Drug analysis in biological fluids

During recent years a consistent number of central nervous system (CNS) drugs have been approved and introduced on the market for the treatment of many psychiatric and neurological disorders, including psychosis, depression, Parkinson disease and epilepsy. Despite the great advancements obtained in the treatment of CNS diseases/disorders, partial response to therapy or treatment failure are frequent, at least in part due to poor compliance, but also genetic variability in the metabolism of psychotropic agents or polypharmacy, which may lead to sub-therapeutic or toxic plasma levels of the drugs, and finally inefficacy of the treatment or adverse/toxic effects. With the aim of improving the treatment, reducing toxic/side effects and patient hospitalisation, Therapeutic Drug Monitoring (TDM) is certainly useful, allowing for a personalisation of the therapy. Reliable analytical methods are required to determine the plasma levels of psychotropic drugs, which are often present at low concentrations (tens or hundreds of nanograms per millilitre). The present PhD Thesis has focused on the development of analytical methods for the determination of CNS drugs in biological fluids, including antidepressants (sertraline and duloxetine), antipsychotics (aripiprazole), antiepileptics (vigabatrin and topiramate) and antiparkinsons (pramipexole). Innovative methods based on liquid chromatography or capillary electrophoresis coupled to diode-array or laser-induced fluorescence detectors have been developed, together with the suitable sample pre-treatment for interference removal and fluorescent labelling in case of LIF detection. All methods have been validated according to official guidelines and applied to the analysis of real samples obtained from patients, resulting suitable for the TDM of psychotropic drugs.

Molecular diversity of Primula apennina and phylogeny of Primula Subsection Euauricula

Primula apennina Widmer is endemic to the North Apennines (Italy). ISSR were used to detect the genetic diversity within and among six populations representative of the species distribution range. High levels of genetic diversity were revealed both at population (PPB = 75.92%, HS = 0.204, Hpop = 0.319) and at species level (PPB = 96.95%, HT = 0.242, Hsp = 0.381). Nei gene diversity statistics (15.7%), Shannon diversity index (16.3%) and AMOVA (14%) detected a moderate level of interpopulation diversity. Principal coordinate and bayesian analyses clustered the populations in three major groups along a geographic gradient. The correlation between genetic and geographic distances was positive (Mantel test, r = 0.232). All together, these analyses revealed a weak but significant spatial genetic structure in P. apennina, with gene flow acting as a homogenizing force that prevents a stronger differentiation of populations. Conservation measures are suggested based on the observed pattern of genetic variability. P. apennina belongs to Primula subsect. Euauricula which includes 15 species distributed on the whole Alps and Apennines. A phylogenetic analysis was carried out using AFLP markers in order both to clarify the relationships among the species of subsection Euauricula that remained unresolved in previous works and to make some hypoteses on their evolutive dynamics. NJ, PCO and BAPS analyses strongly confirmed the monophyly of P. subsect. Euauricula and all the species form strongly supported clades. NJ tree topology suggested a simultaneous fragmentations of ancestral species in a large number of isolated populations that survived in refugia along the unglaciated margins of the Alps in response to the Pleistocene climatic oscillations.

Infezione da virus dell'epatite E nel suino e nell'uomo: aggiornamento sulle metodiche diagnostiche di laboratorio

The first part of the thesis is a brief review on the most important aspects of HEV infection in human and swine, followed by an update on the laboratory techniques currently in use for the diagnosis of HEV infections in humans and animals. The second part refers on the results of two investigations carried out on the presence of HEV infection in swine farms in Toscana and Piemonte and on the presence of HEV infection in pigs and humans in some rural communities in Bolivia. HEV strains isolated from swine herds in Toscana and Piemonte were all included in the genotype 3, showing particular homology with Dutch porcine isolates, Spanish porcine and human isolates and British human isolates. The investigation carried out, with a random sampling, in the province of Cuneo, detected HEV infection with a prevalence of 46% on farms with a number of pigs greater than 500. HEV was detected in pigs and humans in rural communities in Bolivia and all the viral isolate were included in the genotype 3. Aminoacidic homology of human and swine isolates was estimated to be 92%. Results on the development of a Real Time RT-PCR to detect HEV are also reported. The used Real Time RT-PCR protocols, one step and two steps, exhibited good sensitivity to detect several Italian swine HEV strains with high rate of genetic variability.

A candidate gene approach to identify DNA markers associated with meat quality and production traits: investigation of several cathepsin genes in Italian heavy pigs.

The cathepsin enzymes represent an important family of lysosomal proteinases with a broad spectrum of functions in many, if not in all, tissues and cell types. In addition to their primary role during the normal protein turnover, they possess highly specific proteolytic activities, including antigen processing in the immune response and a direct role in the development of obesity and tumours. In pigs, the involvement of cathepsin enzymes in proteolytic processes have important effects during the conversion of muscle to meat, due to their influence on meat texture and sensory characteristics, mainly in seasoned products. Their contribution is fundamental in flavour development of dry-curing hams. However, several authors have demonstrated that high cathepsin activity, in particular of cathepsin B, is correlated to defects of these products, such as an excessive meat softness together with abnormal free tyrosine content, astringent or metallic aftertastes and formation of a white film on the cut surface. Thus, investigation of their genetic variability could be useful to identify DNA markers associated with these dry cured hams parameters, but also with meat quality, production and carcass traits in Italian heavy pigs. Unfortunately, no association has been found between cathepsin markers and meat quality traits so far, in particular with cathepsin B activity, suggesting that other genes, besides these, affect meat quality parameters. Nevertheless, significant associations were observed with several carcass and production traits in pigs. A recent study has demonstrated that different single nucleotide polymorphisms (SNPs) localized in cathepsin D (CTSD), F (CTSF), H and Z genes were highly associated with growth, fat deposition and production traits in an Italian Large White pig population. The aim of this thesis was to confirm some of these results in other pig populations and identify new cathepsin markers in order to evaluate their effects on cathepsin activity and other production traits. Furthermore, starting from the data obtained in previous studies on CTSD gene, we also analyzed the known polymorphism located in the insulin-like growth factor 2 gene (IGF2 intron3-g.3072G>A). This marker is considered the causative mutation for the quantitative trait loci (QTL) affecting muscle mass and fat deposition in pigs. Since IGF2 maps very close to CTSD on porcine chromosome (SSC) 2, we wanted to clarify if the effects of the CTSD marker were due to linkage disequilibrium with the IGF2 intron3-g.3072G>A mutation or not. In the first chapter, we reported the results from these two SSC2 gene markers. First of all, we evaluated the effects of the IGF2 intron3-g.3072G>A polymorphism in the Italian Large White breed, for which no previous studies have analysed this marker. Highly significant associations were identified with all estimated breeding values for production and carcass traits (P<0.00001), while no effects were observed for meat quality traits. Instead, the IGF2 intron3-g.3072G>A mutation did not show any associations with the analyzed traits in the Italian Duroc pigs, probably due to the low level of variability at this polymorphic site for this breed. In the same Duroc pig population, significant associations were obtained for the CTSD marker for all production and carcass traits (P < 0.001), after excluding possible confounding effects of the IGF2 mutation. The effects of the CTSD g.70G>A polymorphism were also confirmed in a group of Italian Large White pigs homozygous for the IGF2 intron3-g.3072G allele G (IGF2 intron3-g.3072GG) and by haplotype analysis between the markers of the two considered genes. Taken together, all these data indicated that the IGF2 intron3-g.3072G>A mutation is not the only polymorphism affecting fatness and muscle deposition in pigs. In the second chapter, we reported the analysis of two new SNPs identified in cathepsin L (CTSL) and cathepsin S (CTSS) genes and the association results with meat quality parameters (including cathepsin B activity) and several production traits in an Italian Large White pig population. Allele frequencies of these two markers were evaluated in 7 different pig breeds. Furthermore, we mapped using a radiation hybrid panel the CTSS gene on SSC4. Association studies with several production traits, carried out in 268 Italian Large White pigs, indicated positive effects of the CTSL polymorphism on average daily gain, weight of lean cuts and backfat thickness (P<0.05). The results for these latter traits were also confirmed using a selective genotype approach in other Italian Large White pigs (P<0.01). In the 268 pig group, the CTSS polymorphism was associated with feed:gain ratio and average daily gain (P<0.05). Instead, no association was observed between the analysed markers and meat quality parameters. Finally, we wanted to verify if the positive results obtained for the cathepsin L and S markers and for other previous identified SNPs (cathepsin F, cathepsin Z and their inhibitor cystatin B) were confirmed in the Italian Duroc pig breed (third chapter). We analysed them in two groups of Duroc pigs: the first group was made of 218 performance-tested pigs not selected by any phenotypic criteria, the second group was made of 100 Italian Duroc pigs extreme and divergent for visible intermuscular fat trait. In the first group, the CTSL polymorphism was associated with weight of lean cuts (P<0.05), while suggestive associations were obtained for average daily gain and backfat thickness (P<0.10). Allele frequencies of the CTSL gene marker also differed positively among the visible intermuscular extreme tails. Instead, no positive effects were observed for the other DNA markers on the analysed traits. In conclusion, in agreement with the present data and for the biological role of these enzymes, the porcine CTSD and CTSL markers: a) may have a direct effect in the biological mechanisms involved in determining fat and lean meat content in pigs, or b) these markers could be very close to the putative functional mutation(s) present in other genes. These findings have important practical applications, in particular the CTSD and CTSL mutations could be applied in a marker assisted selection (MAS) both in the Italian Large White and Italian Duroc breeds. Marker assisted selection could also increase in efficiency by adding information from the cathepsin S genotype, but only in the Italian Large White breed.

L’allevamento della sogliola comune (Solea solea L.): riproduzione controllata e nuovi sviluppi di nutrizione e alimentazione in fase larvale e giovanile

In questa tesi sono illustrate alcune sperimentazioni finalizzate alla standardizzazione del ciclo produttivo della sogliola comune (Solea solea) in cattività. E’ stato creato un parco di riproduttori selvatici ed è stata standardizzata la riproduzione ad un livello compatibile con la realtà produttiva del settore. Indagini genetiche di assegnazione parentale hanno evidenziato come alcuni esemplari siano stati predominanti negli accoppiamenti e nel conseguente contributo alla generazione della prole. Ciò ha determinato una diminuzione della variabilità genetica dei discendenti. La composizione quali-quantitativa degli acidi grassi delle uova è stata correlata con la sopravvivenza larvale nel corso di un’intera stagione riproduttiva. Tale composizione non ha subito importanti variazioni su scala temporale e sembra essere stata influenzata dall’alimentazione somministrata ai riproduttori nel periodo precedente alla riproduzione. Le analisi di interazione tra momento riproduttivo e qualità delle uova hanno confermato che è stato possibile ottenere uova di buona qualità in termini di sopravvivenza larvale nel corso di tutta la stagione riproduttiva. Larve di sogliola sono state svezzate precocemente 13 giorni dopo la schiusa riducendo l’impiego di cibo vivo a favore di micro diete commerciali. Tale svezzamento ha ridotto le performance di accrescimento, ma non la sopravvivenza e lo sviluppo della metamorfosi quando comparati ad un trattamento standard. La riduzione del cibo vivo ha ottimizzato i costi di produzione e migliorato l’igiene in vasca. L’ontogenesi di precursori di enzimi digestivi è stata determinata tramite PCR quantitativa. I risultati di espressione di tripsinogeno, chimotripsinogeno e amilasi hanno mostrato come tali enzimi rivestano un ruolo chiave nei processi digestivi delle prime fasi larvali. Esemplari giovanili hanno ottenuto un significativo maggiore indice di accrescimento e migliore indice di conversione quando alimentati con diete sperimentali contenenti un elevato tenore proteico. Un aumento dell’incidenza di vacuoli lipidici a livello epatico è stato osservato all’aumentare del tenore proteico della dieta.

Taxonomy, phylogeny and reproductive ecology of Gentiana lutea L.

This research focuses on taxonomy, phylogeny and reproductive ecology of Gentiana lutea. L.. Taxonomic analysis is a critical step in botanical studies, as it is necessary to recognize taxonomical unit. Herbarium specimens were observed to assess the reliability of several subspecies-diagnostic characters. The analysis of G. lutea genetic variability and the comparison with that of the other species of sect. Gentiana were performed to elucidate phylogenetic relationships among G. lutea subspecies and to propose a phylogenetic hypothesis for the evolution and the colonization dynamics of the section. Appropriate scientific information is critical for the assessment of species conservation status and for effective management plans. I carried out field work on five natural populations and performed laboratory analyses on specific critical aspects, with special regard to G. lutea breeding system and type and efficiency of plant-pollinator system. Bracts length is a reliable character to identify subsp. vardjanii, however it is not exclusive, hence to clearly identify subsp. vardjanii, other traits have to be considered. The phylogenetic hypotheses obtained from nuclear and chloroplast data are not congruent. Nuclear markers show a monophyly of sect. Gentiana, a strongly species identity of G. lutea and clear genetic identity of subsp. vardjanii. The little information emerging from plastid markers indicate a weak signal of hybridization and incomplete sorting of ancestral lineages. G. lutea shows a striking variation in intra-floral dichogamy probably evolved to reduce pollen-stigma interference. Although the species is partially self-compatible, pollen vectors are necessary for a successful reproduction, and moreover it shows a strong inbreeding depression. G. lutea is a generalist species: within its spectrum of visitors is possible to recognize "nectar thieves" and pollinators with sedentary or dynamic behaviour. Pollen limitation is frequent and it could be mainly explained by poor pollen quality.

Isolamento, variabilità ed espressione del gene PIMT in Girasole

Il vigore è un aspetto rilevante della qualità delle sementi, strettamente connesso al loro status fisiologico, al genotipo e alle condizioni di stoccaggio. Stress ossidativi e danni alle macromolecole sono alla base del deterioramento del seme, che è equipaggiato con sistemi protettivi e di riparazione. Uno di questi coinvolge l’L-isoaspartil metiltransferasi (PIMT) che ripara i misfolding proteici catalizzando la riconversione in aspartato dell’isoaspartile anomalo accumulato. Scopo di questo studio era valutare il possibile ruolo del meccanismo di riparazione di PIMT nel vigore del seme in girasole. Per questo il relativo gene è stato isolato e caratterizzato, la variabilità allelica determinata su un campione di linee inbred e l’espressione genica misurata in risposta all’invecchiamento accelerato (aging) e al priming. La sequenza codificante ottenuta è costituita da 4 esoni e contiene i 5 domini caratteristici delle metiltransferasi. Il gene mostra elevata similarità con gli ortologhi vegetali e scarsa diversità nucleotidica nei genotipi coltivati rappresentativi della variabilità della specie. Nella sequenza aminoacidica, comunque, sono state rinvenute tre sostituzioni che potrebbero influenzare la funzionalità enzimatica. Dal punto di vista fisiologico i genotipi considerati hanno esibito notevole variabilità di risposte ai trattamenti, sia in termini di vigore che di espressione genica. Aging e priming hanno prodotto generalmente gli effetti attesi, rispettivamente negativi e positivi, sulla germinabilità e sulla sua velocità. In generale l’espressione di PIMT è risultata massima nel seme secco, come riportato altrove, e ridotta dall’aging. Anche il priming ha diminuito l’espressione rispetto al seme quiescente, mentre il suo effetto dopo l’aging è risultato genotipo-dipendente. Tuttavia, nelle condizioni descritte, non si sono evidenziate correlazioni significative tra vigore ed espressione di PIMT, tali da suggerire un chiaro ruolo di questo meccanismo nella qualità fisiologica del seme.

Conservation genetics of European wildcat (Felis silvestris silvestris): a wide and integrating analysis protocol for admixture inferences and population structure

Introgression of domestic cat genes into European wildcat (Felis silvestris silvestris) populations and reduction of wildcats’ range in Europe, leaded by habitat loss and fragmentation, are considered two of the main conservation problems for this endangered feline. This thesis addressed the questions related with the artificial hybridization and populations’ fragmentation, using a conservation genetics perspective. We combined the use of highly polymorphic loci, Bayesian statistical inferences and landscape analyses tools to investigate the origin of the geographic-genetic substructure of European wildcats (Felis silvestris silvestris) in Italy and Europe. The genetic variability of microsatellites evidenced that European wildcat populations currently distributed in Italy differentiated in, and expanded from two distinct glacial refuges during the Last Glacial Maximum. The genetic and geographic substructure detected between the eastern and western sides of the Apennine ridge, resulted by adaptation to specific ecological conditions of the Mediterranean habitats. European wildcat populations in Europe are strongly structured into 5 geographic-genetic macro clusters corresponding to: the Italian peninsular & Sicily; Balkans & north-eastern Italy; Germany eastern; central Europe; and Iberian Peninsula. Central European population might have differentiated in the extra-Mediterranean Würm ice age refuge areas (Northern Alps, Carpathians, and the Bulgarian mountain systems), while the divergence among and within the southern European populations might have resulted by the Pleistocene bio geographical framework of Europe, with three southern refugia localized in the Balkans, Italian Peninsula and Iberia Peninsula. We further combined the use of most informative autosomal SNPs with uniparental markers (mtDNA and Y-linked) for accurately detecting parental genotypes and levels of introgressive hybridization between European wild and domestic cats. A total of 11 hybrids were identified. The presence of domestic mitochondrial haplotypes shared with some wild individuals led us to hypnotize the possibility that ancient introgressive events might have occurred and that further investigation should be recommended.